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Abstract. By virtue of numerical arguments we study a bifurcation phenomenon occurring for a class of 
minimization problems associated with the quasi-linear Schrodinger equation. 



1. Introduction 

Various physical situations are described by quasi-linear equations of the form 



(1.1) 



i0t + + 0A|(^|2 + = in (0, oo) x R^, 



(j){0,x) = (f>o{x) 



where 1 < p < 11, i stands for the imaginary unit and the unknown </> : (0, oo) x — C is a complex 
valued function. For example, it is used in plasma physics and fluid mechanics, in the theory of Heisenberg 
ferromagnets and magnons and in condensed matter theory. See e.g. the bibliography of [6]. Motivated by 
the classical stability results of the semi-linear Schrodinger equation 



(1.2) 



i(l)t + A(t)+\<j)\P-^(t)^0 in(0,oo)xR^ 
0(0, x) = 00 (a;) inM^ 



namely the stability of ground states (least energy solutions of —Am + uju = \u\p~^u, w > 0) for 1 < p < | 
and their instability for p > | (see [2,5]), an interesting and physically relevant question for equation (1.1) 
is the orbital stability of ground state solutions of 

(1.3) - Au - uAu'^ + oju = \u\P-^u inR^. 

When 1 < p < it is conjectured in [6] that the ground states are orbitally stable. However, in [6] this 
result was not proved. Instead, it was considered the stability issue for the minimizers of the problem 

(1.4) ^(c) = inf ^(u), 

ll"ll_f,2(p_3)=C 

where 1 < p < ^ and the energy functional defined on X = {ue H\m.^) : u\Du\ e L^{R^)} by 



- f {l + 2u^)\Du\^dx — f \u\P+^dx. 

2 Jk3 p+l Jus 



(1.5) ^(u) 2 

This problem, which looks interesting by itself, can be seen as a useful tool for a first attempt towards the 
understanding of orbital stability of ground states of (1.3) for fixed lu > 0. Denoting by Q{c) the set of 
solutions to (1.4), in [6] the authors prove that if 1 < p < ^ and c > is such that m(c) < 0, then Q{c) is 
orbitally stable (see [6] for the definition). Concerning (1.4), we learn [6] that the following facts hold: 

Proposition 1.1 (CJS, [6]). The following properties hold. 

(1) Assume that 1 < p < Then ^(c) > — oo for every c > 0. 

(2) Assume that p = Then ^(c) > — oo for c > small and ^(c) = — oo for c > large. 

(3) Assume that ^ < p < 11. Then ^[c] = — oo for every c > 0. 
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(4) Assume that 1 < p < ^ and let c > be such that ^{c) < 0. Then problem (1.4) admits a positive 
minimizer Uc €z X which is radially symmetric and radially decreasing. Moreover any solution 
Uc G X of (1.4) satisfies 

(1.6) - AUc - UcAul + XcUc = \Ucf~^Uc 

for some Ac > 0. 

(5) Assume that 1 < p < Then ^{c) < for every c > 0. 

(6) Assume that | < p < Then ./#(c) < for every c > 0. 

(7) Assume that | < p < Then there exists 

C(j = c(p) > 

such that: 

(a) If c < Cf^ then ^(c) = and ^{c) does not admit a minimizer. 

(b) If c > c^ then ^(c) < and ^(c) admits a minimizer. 

The motivation for the formulation of the following problems is mainly related to point (7) in the 
statement of Proposition 1.1. Notice that the value of C(j(p) can be characterized as follows 

cj(p) = inf{c > : .Jf{c) < 0}, 

This could help while trying to numerically compute the value of cu (p) . The bifurcation value cj (p) which 
appears in the previous Proposition 1.1 is obtained in [6] by an indirect argument by contradiction and 
thus it is not explicitly available for calculation through a given formula. Hence, on these basis, it seems 
natural to formulate the following problems: 

Problem 1.2. Provide some lower and upper bounds of ^{c) for c > and p e ^). 
Problem 1.3. Numerically compute or provide bounds for the map cj : [|, ^) — > (0, +oo). 
Problem 1.4. Numerically compute the solutions to .^(c) for c > cj(p) with p G ^). 
For the corresponding, more classical [5], semi-linear minimization problem 

^si(c) = inf ^si{u), S'4u) = i / \Du\''dx - f \u\P+^dx, K p < ^, 

ll«lli2(E3)=C 

there is no bifurcation phenomena, namely C|j(p) = for every 1 < p < ^ and Problem 1.4 was studied 
in [4] by arguing on a suitable associated parabolic problem in order to decrease initial energies computed 
on Gaussian initial guesses. An important point both for analytical and numerical purposes is the fact that 
minimizers of ^si(c) have fixed sign, are radially symmetric, decreasing and unique, up to translations and 
multiplications by ±1. In principle, the uniqueness is used in [4] to justify that the numerical algorithm 
really provides the solution to the minimization problem ^si(c). To show the uniqueness of solutions to 
^si(c) one can argue as follows. By the result of [9], for any A > there exists a unique (up to translations) 
positive and radially symmetric solution r = rx '■ M'^ — !• R of 

-Ar + Xr = rP in R^ 

In turn, given Ai, A2 > 0, if ri, r2 : K'^ M denote, respectively, positive radial solutions of 

(1.7) -Ari + Airi=rP in M^ -Ar2 + A2r2 = in 
there exists some point ^ G M"^ such that 

(1.8) r2(x + = ^(72;), ^---{^y^ ^'-^i^)^- 

Let now ri and r2 be two given solutions to the minimization problem ^si(c) and let — Ai and — A2 be the 
corresponding Lagrange multipliers. By virtue of [5, Theorem II. 1, ii)], Ai, A2 > and ri, r2 are solutions 
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of (1.7), are radially symmetric, radially deereasing and with fixed sign. In turn, up to multiplication by 
±1, ri,r2 > and by (1.8) it holds 



C= il''2ilL2(R3) = lk2(-+0llL2(R3) 7 J^Jiix)dx C^i J , 

yielding in turn n'^j"^ — 1. By the definition of 7 and /i in (1.8), we get Ai = A2 and 7 = /i = 1, yielding 
from (1.8) as desired r\ = r2, up to a translation. 

On the contrary, it is not currently known that minimizers r £ X of the quasi-linear minimization 
problem (1.4) are unique, up to translations and multiplications by ±1, although it is conjectured that this 
is the case. If r £ X is a given minimizer for (1.4), arguing as in [6] it is possible to prove that it has fixed 
sign, so that, up to multiplication by —1, we may assume r > 0. Then, r is radially symmetric and radially 
decreasing, see for instance the main result of [12]. Given ri and r2 two solutions to the minimization 
problem ^(c) we have that — Ai and — A2 are the corresponding Lagrange multipliers and Ai,A2 > in 
light of [6, Lemma 4.6] 

(1.9) - Ari - r,l^r1 + Xir, = rf in M^ i = l,2. 

With respect to the semi-linear case, the main problem is the identification of Lagrange multipliers, which 
cannot be inferred as in the semi-linear case. In fact, let us consider as above the rescaling for ri by 



w{x) = ^lrl{-1x), ^■.^(^^\ ^.-(^^ 



Then this yields ||w|1^2(r3) = (A2/A1) ''^('p-i'' c and 



-Aw - (Ai/A2)''-i wAw^ -I- A2W = inM^. 

Although there are recent uniqueness results for the positive radial solutions to (1.9), due to the presence 
of the residual coefficient (Ai/A2)^/'-^~^-', we cannot infer as before that w{-) = r2(- + ^) for some point 
^ S R^. In the course of the next section, we shall compute the ground state under a conjectured property 
(indeed true in the semi-linear case discussed above) stated in the following 



Conjecture 1.5. Assume that ri,r2 > are radial decreasing solutions to (1.9) with Ai, A2 > 0, £{ri) < 
and ||7'i||^2(R3) = c for i = 1,2. Then Ai = A2 and r2 = 7'i(- l^r some ^ G R'^. 

Consequently, given c > 0, assume that there exist r : R'^ — > R and A > such that 

r = pda;]) > 0, p' < 0, lklli2(R3) = c, f (r) < 0, -Ar - rAr^ + Ar = in R^ 

We know that this happens to be the case under the assumption that ^(c) < 0. Then r is the unique 
solution to problem ^(c), up to translations and multiplication by ±1. 
In figure 1, we compared the shape of the solutions to 

(1.10) ^^(70) ~ inf ^i;(u), S^{u):=l-l {I + 2§u^)\Du\^ dx - \ f \u\^dx. 

"S-^ 2 Jr3 3 Jr3 

Il"lll2,,:3,=70 

with '8 = 1 (quasi-linear case) and ■(? = (semi- linear case). Roughly speaking, the term — mAu^ produces 
an additional diffusive contribution which tends to squeeze the bump down against source effects. 

2. Results 

Concerning Problem 1.2, we have the following (see also Fig. 2) 

Proposition 2.1. The following properties hold. 
(1) For every ^ < p < ^ there holds 

O O 2 r, 2+1 3 3(P+1) 

280 5 3C^ 3(p-l) 2V2C 2 7r2 ' 



VoO: ^{c)< M {a^—+a''—^—-a 



'^>o\ 4 8V27r3/2 {p+l)i 
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Figure 1. Comparison between the solution (a section of the absolute value squared) of (1.10) 
with i9 = 1 (quasi- linear case) and with i9 = (semi-linear case). The additional diffusive term 
present in the quasi-linear case tends to produce squeezing effects. 



(2) For every I <p < ^, if 



/3^rV3(3/2)' 

being Ssoh the best Sobolev constant for the embedding of H^{S?) into -L'(M^), there holds 

(3) For p = |, setting A := ^^j^, B := and cj := ( 2^233/2 )^^'^, there holds 

/n particular, ^{c)\^^t_ < — Mc^'^/^, /or ewer?/ c > large and some M > 0. 

(4) If p= setting 

c, (16/3/C 13)3/2 ^ 73^ 452^^^ 

we /laue ^{c) ~ /or ewery c < C[, and the infimum ^(c) is not attained. 

(5) If p — setting 

C . ^^3^2 ~ 00. U», 

we /lave ^(c) = —00 /or every c> c^ . 

Proof. Properties (1), (3) and (5) easily follow by the arguments in Section 3 and direct computations. 
Properties (2) and (4) need bounds from below and can be justified as follows. The best Sobolev constant 
S'sob is computed through the formula contained in [13]. Concerning (4), by Holder and Sobolev inequalities, 
for It € X we have 

(2.1) / \u\^'^/^dx<( [ / lul^^dx)'^"^ < /Ci3c2/3( / \u\^\Du\^dx 
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where we have used the fact that (2* = 6 is the critical Sobolcv exponent in 

/ \u\^^dx= [ {u^f'dx, [ \D{u^)\^dx = 4 [ \u\^\Du\^dx. 
From inequahty (2.1) we infer 

r 3^13 , r 

£{u)> \u\''\Du\^dx-—^^'^ luHz^upd^, 

JR3 16 Jr3 



which yields £{u) > for every u € X and any c < and hence, in turn, the desired conclusion. In a 
similar fashion, concerning (2), if p < 13/3, using Holder and Sobolev inequalities for any u ^ X we have 



3p-3 



(2.2) / \u\P+^dx<ICpC^( [ \u\'^\Vu\^dx) 

yielding immediately that 

£{u)> I \u\^\Vu\^dx-^c^{ I \u\^\Vu\^dx\^ . 

Jr3 p+l \ ^ 

Since p < 3/13 it follows that < i. in turn, the function cj : [0, +oo) ^ M 

. . /Cn 11-P 3p-3 

Uj(s) = S —C 10 S 10 , s > 0, 

^ ' p+l 

always admits a (negative) absolute minimum point at a point <p,c > which can be easily computed, 
yielding the desired assertion by the arbitrariness of u e X. □ 

Concerning Problem 1.3 we shall provide an upper bounding profile for the values of c\j{p) and give indi- 
cations showing that very likely C|)(p) remains in a small lower neighborhood of this profile. As p increases 
from 7/3 up to a certain value po ^ 3.3 the bounding profile is increasing, reaching values around c = 250. 
Then, after po it decreases. 

Concerning Problem 1.4, under the conjectured uniqueness result we compute the ground state solutions 
for some values of c greater than the upper bounding profile for the values of c^{p)- In the case c < c^{p), 
roughly speaking, if (u„) C X is an arbitrary minimizing sequence for problem ^(c), since we know that 
^(c) is not attained, (u„) cannot be strongly convergent, up to a subsequence, in L^(R'^) and in turn in 
{JS? ) . Then, by virtue of Lions's compactness-concentration principle, only vanishing or dichotomy might 
occur, in the language of [10,11]. On the other hand, it was proved in [6, Theorem 1.11] that dichotomy can 
always be ruled out. In conclusion, the only possibility remaining for a minimizing sequence is vanishing, 
precisely: 

for every i? > 0: lim sup / v^dx = 0, 

"■^°°yeR3 JB(y,R) 

where -B(t/, R) denotes the ball in 'M? of center y and radius R. In particular, fixed any bounded domain C in 
R"^, imagined for instance as the computational domain, the sequence (it„) cannot be essentially supported 
into C, being for i?o = diam(C) 

/ u\dx > c — sup / u^dx > c — sup / u^dx > c/2, 

JR3\C yecJB(y,Ra) yeM.^ J B(y,Ro) 

for n sufficiently large. This means that (u„) (in particular any numerically approximating ^(c)) tends 
to spread out of any fixed (computational) domain C. For instance, the sequence {gc,aj ) with aj — > as 
j — oo, where gc,a is defined in Section 3 is such that ||(7c,o-j 11^2(^3) = c and £{gc,aj) — as j — > oo and 
hence, for c < C[i(p), being ^(c) = 0, it is a minimization sequence. It shows vanishes, since 



for all i? > and j > 1. 



sup / g-'^^ dx < CR'^a'^, 

'Buiv)' ' ' 
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Figure 2. Upper and lower bounds of .^(c) according to the estimates obtained Proposition 2.1 
in the particular cases p — 7/3 (sharp bounds) and p — 10/3 (less sharp bounds). 



3. Numerical approximation 

Instead of a direct minimization of the energy fmictional (1.5) (see, for instance, [3]), as seen in [1,4], it 
is also possible to find a solution of (1.6) by solving the parabolic problem 



(3.1) 



dtu ^ Au + uAu'^ + \u\P ^u~X{u)u 
u(0) = uo, lko|ji2(R3) = c 



up to the steady-state Uc, where X{u) is defined by 



X{u) = 



/r 



This approach in known as continuous normalized gradient flow. It is easy to show that the energy associated 
to the solution u(t) decreases in time, whereas the L^-norm is constant. 
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Figure 3. Values of the infimum (Gaussian) energy with respect to c (left) and values of Cg (right): 
from p = 7/3 to p = 12/3 they are 69.510, 150.65, 212.91, 241.01, 225.42, 170.33, respectively. 



In order to find a "good" initial solution uq, that is a function with a negative energy, we can consider 
the family of Gaussian radial functions 



9c.air) 



I15c,(t|1l2(r3) — c, 



a>0 



and minimize the energy S'{gc^a-), which, for a given p and c, can be computed analytically as 



(3.2) 



, 3c 



3c2 



2V2, 



P + 1 3 3(p + l) 
C 2 7r2 4 



8\/27r3/2 



(p+1) 2 



with respect to the parameter a. If the infimum value for the energy is zero, we increase the value of c 
and look again for the infimum energy. We proceed with increasing values of c until we find a Cg such that 
the minimum energy with respect to a, corresponding to a value a, is negative. This is possible, since we 
consider a discrete sequence of increasing values for c in order to test the negativity of the energy. Such a 
Cg is clearly an upper bound for the desired value cj. For the range p = 7/3, 8/3, 9/3, 10/3, 11/3, 12/3 we 
obtain the values reported in Figure 3 (right). 
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Now, we arc ready to look for values c < Cg, for which the steady-state of (3.1) has a negative energy. 
To this purpose, we choose uq = gc,a, where a = a ■ c/cg. The meaning of this choice is the following: since 
c < Cg; the infimum value of the energy attained by a Gaussian function is zero and corresponds to the 
limit case a ^ 0. We instead select < a < a with the idea that gc^a is a good initial value, because close 
to Qcg.a, that is the optimal element in the Gaussian family. With this choice, clearly we have 'S'{uq) > 0. 

In order to fix once and for all the computational domain in such a way that it does not depend on a, 
we scale the space variables by a and the unknown u in order to have unitary norm, that is 



We end up with 

(3.3) 

where 



V ccr^Vc{t, a-) ~ u{t, •) 



/k3 ""c 



1 f „JiiA , o„-3„2N|n„ |2 _ c(eo-^) ^ ; I ip+i 



The corresponding energy is 

S(u) = E{v,) = ^ / ca'il + 2ca''vi)\Dv^, 

2 Jr3 P + 1 

We solve equation (3.3) up to a final time T for which 

\\a^Avc + ca'^VcAvl + {ca^)^\vc\P-^Vc + vivc)vc\\mM^) < tol 

where tol is a prescribed tolerance to detect the approximated steady-state. As already done in [4], we 
apply the exponential Runge-Kutta method of order two (see [8]) to the spectral Fourier decomposition in 
space of (3.3). The embedded exponential Eulcr method gives the possibility to derive a variable stepsize 
integrator, which is particularly useful when approaching the steady-state solution, allowing the time steps 
to become larger and larger. For our numerical experiments, we used the computational domain [—5,5]"^, 
the regular grid of 64"^ points, a tolerance for the local error (in the norm) equal to 10~® and the 
steady-state detection tolerance tol = 10~^. 

The solution Vc(T) is then considered an approximated steady-state solution. If its energy is negative 
and it is radially symmetric and decreasing, then wc conclude it is a minimum for ^(c). Therefore, c > C(| 
is the current upper bound for cj (blue circle in Figure 4). 

On the other hand, if c < cu, then ^(c) — and the infimum is approximated by flatter and flatter 
functions. This is perfectly clear in the case one restricts the search among Gaussian functions, for which 
<^(ffc,cr) for (T — > 0. This situation can be recognized in the numerical experiments because the essential 
support of Vc{t) during time evolution tends to grow and to spread out the computational domain (green 
plus in Figure 4) , which was chosen in such a way to comfortably contain the essential support of the initial 
solution Vc{0). In this case, we apply a bisection algorithm on the values c and Cg (the Gaussian upper 
bound for cj) in order to flnd a tighter upper bound for cj, since c < cj and Cg > cj. 

In the numerical experiments wc encountered another situation: starting from an initial solution with 
positive energy, it was possible to find a radially symmetric and decreasing approximated steady-state 
solution, whose support was perfectly contained into the computational domain and with a positive energy 
(red star in Figure 4) . This solution is not a solution with minimum energy, because for a small enough it 
is possible to find a Gaussian function gc^a with smaller energy. In our numerical experiments, we observed 
this behaviour, for a given p and with the tolerances described above, for the values of c between the first 
value for which the essential support of the solution spread out the computational domain and the current 
upper bound for cj. We notice that we were not able to obtain solutions with positive energy in the limit 
case p = 7/3. 

Overall, from the numerical experiments we found out that the higher is p the more difficult is to find a 
value of c < Cg for which the steady-state has a negative energy and can thus be considered a solution of 
minimum energy. 
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100 - 



Figure 4. Values of Cg (Gaussian upper bound of cj), Cm (values for which E{vc„^{T)) < 0) and Cp 
(values for which E{vc-p (T)) > 0). For instance, for the case p = 9/3 and c = 0.9875-Cg ~ 210.25 we 
found an approximated steady-state with energy E{vc{T)) f» —1.30 ■ 10^^ (blue circle), whereas 
with c = 0.975 • Cg ~ 207.59, we found an approximated steady-state with energy E{vc{T)) ~ 
6.97 • 10"^ (red star). 
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